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MODEL FOR TESTING IMMUNOGENICITY OF PEPTIDES 

Governtnen t Interest 

The invention described herein may be manufactured, 
licensed and used by or for governmental purposes without the 
payment of any royalties to us thereon. 

Cross Reference 

This application is a continuation-in-part of U.S. Patent 
application Serial No. 08/064,559, filed May 21, 1993, and the 
present application incorporates U.S. Patent Application 
Serial No. 08/064,559 in its entirety by reference. 

Field of the Invention: 

This invention relates to a means of predicting potential 
of a peptide for eliciting immune response. 

Background of the Invention: 

Among the numerous steps required for an immunological 
response to occur is the presentation of the antigen by 
macrophages to the B-cell or T-cell. This presentation is 
mediated by the Class I and Class II major histocompatibility 
complex (MHC) molecules on the surface of 1 the cell. The MHC 
molecules hold antigens in the form of the peptide fragments 
and together with the receptor molecule on the T-cells, form a 
macromolecular complex that induces a response in the T-cell. 
Therefore, a necessary step in an immune response is the 
binding of the antigen to the MHC. 
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Recent single crystal X-ray structures of human and 
murine Class I MHC's have been reported. Analysis of these 
crystal structures have shown that antigenic peptides lie in 
the so-called binding cleft for presentation to the T-cell. 
This cleft is formed by a l and a 2 domains and by ]8- strands from 
each domain forming the floor. Furthermore, the sequence 
polymorphism among Class I molecules can result in alterations 
of the surface of the cleft forming different pockets. 
Peptide side chains may insert into these pockets. Thus, 
different pockets may interact with different side chains. 
This implies the mechanism for the peptide specificity of 
Class I MHC's. Peptides . bound to the Class I MHC's in the 
crystal structures were found to have both the amino and 
carboxy termini tightly held by the MHC. There were few 
interactions near the middle of the cleft. Hence the bound 
peptide is allowed to bend slightly in the center. The 
observed binding mode helped to explain the apparent partial 
specificity of peptide sequence and the allowed variation in 
peptide length found among peptides isolated from Class I 
MHC' s . 

The precise mode of binding of peptides to Class II MHC 
molecules is less clear. While a single crystal X-ray 
diffraction structure for the HIA-DR1 MHC has been shown, the 
coordinates have remained unavailable. However, currently 
available theoretical and experimental results help form a 
hypothesis that the binding of a peptide to Class II MHC is 
similar to that observed with Class I. First, it is noted 
that the Class II binding cleft is structurally similar to 
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that of Class I. This was concluded based upon a sequence 
analysis of 26 Class I and 54 Class II amino acid sequences. 

Unlike with Class I molecules, self -peptides isolated 
from murine I-A b and I-E b , from murine I-A d and from human 
5 HLA-DR1 molecules were found to be varied in size (13 to 25 

residues long) . The peptides isolated from the murine I-A b and 
I-E b molecules had heterogenous carboxy termini while those 
from I-A d and HLA-DR1 had ragged termini at both ends. The 
varying lengths indicate that the amino and carboxy termini of 
10 the peptides were not critical for the binding. One or both 
termini may protrude from the binding site and be available 
for further processing. The residues critical for binding 
were proposed to be at the ends of the peptide as opposed to 
the center. 

15 

Summary of the Invention: 

It is the purpose of this invention to provide a method for 
preliminary screening of peptides for ability to elicit an immune 
response . Structural homology techniques were used to model a 
20 receptor (the Class II MHC is exemplified) . This model makes it 

possible to preliminarily screen peptides for antigenic properties. 
By modifying the peptide to "fit" into the receptor it is possible 
to identify methods of rendering non- immunogenic peptides 
immunogenic . 

25 The preliminary screening of peptides for immunogenic ity 

comprises the steps of (1) creating a molecular model of a receptor 
followed by minimizing the model created, 2) modeling a peptide to 
be tested and minimizing the model of the peptide, then testing the 
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fit of the model of the peptide into the model of the receptor to 
produce a composite minimized receptor /minimi zed peptide model - 
Upon finding an acceptable fit, the peptide may then be screened by 
a binding assay for actual binding to Class II MHC as a further test 
for immunogenicity . 

It has been found that when the model of the peptide can not be 
fitted into the model of the receptor, the peptide will lack 
immunogenicity. While not all peptide models which can be made to 
"fit" into to model of the receptor will be effective as immunogens, 
the screening methods of the invention may make it possible to avoid 
undue biological testing of inappropriate peptides. By using the 
model, it is also possible to alter peptides to accommodate the 
receptor. Hence, the invention has both predictive and drug design 
applications . 

Brief Description of the Figures: 

Fig. 1 shows the HliA-aw68 a, and ot 2 domains with DR1 oc l and 0 X 

domains . 

Figs. 2-30 are a printout of the minimized coordinates of the 
receptor . 

Figs. 31 and 32 shows the effects of various peptides 
inhibiting the binding of labeled hemagglutinin in a competitive 
binding assay. 

Detailed Description of the Inven tion: 

In order to understand and better predict peptide interaction 
with Class II MHC's and as an aid for synthetic peptide vaccine 
design, a structural homology model of HLA-DRl molecule was made 
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using the Class I HLA-aw68 as a reference molecule. For purposes of 
this analysis, numerous conserved residues were aligned leading to a 
proposed three-dimensional model for the Class II structure very 
similar to that of Class I . This model retained the overall 
5 conformation of a Class I MHC and agreed with a considerable amount 
of the published data. Furthermore, peptides shown to bind to DR1 
were docked in the binding cleft of the model and analyzed. The 
results agree with the experimental binding data presented here. 
Hence, it is shown that the structural homology model reported here 

10 is useful for screening Class II MHC functionality. 

It had been hypothesized that few peptide residues may be 
required for binding to DR1 . By substituting residues into the ■ 
influenza hemagglutinin 307-319 T-cell epitope (HA) it had been 
determined that a single tyrosine at 3 08 was required for binding. 

15 A synthetic peptide with the tyrosine at position 308 and a lysine 
at 315 was found to bind DR1 as well as the native peptide. Hence, 
it was concluded that few peptide residues determine the high 
affinity binding to DR1 . 

The peptides produced according to the present invention may be 

20 used alone or chemically bound to another peptide and/or carrier in 
order to elicit an immune response. An immune response is elicited 
by administering a peptide to an animal in an effective dose and by 
an effective route of administration. Typically the peptide will be 
administered with an immunologically acceptable carrier. The routes 

25 of administration, dosages, times between multiple administrations 
will be based on the particular peptide and are standard operations 
of those skilled in the art. 
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Of particular interest are peptides from pathogenic 
microorganisms and neoplasms. In such an example, a vaccine may be 
formed with the peptide and any known immunological carrier and may 
be administered prophylactically or therapeutically. The immune 
5 response may be elicited for a number of reasons other than for 

prophylaxis or therapy such as increasing antibody production for 
the harvesting of antibodies, or increasing specific B-cell or T- 
cell concentration for the production of hybridomas or cellular 
therapy . 

10 The choice of host animals is limited only to those capable of 

an immune response. Preferred hosts are mammals, more preferred are 
humans . 

The vaccine may contain plural peptides with each peptide 
corresponding to the same or different antigens. The peptides may 

15 be used unbound or they may be chemically bound to another peptide 
or an unrelated protein or other molecule . A preferred vaccine 
preparation contains a plurality of peptides chemically bound to a 
larger more immunogenic peptide. 

The peptide may be adsorbed, bound or encapsulated in a 

20 biodegradeable microsphere, microcapsule, larger carrier or a 

combination of these . The carrier may have a slow or controlled 
release property thereby releasing the peptide under appropriate 
conditions and times for enhanced immunization. This is 
particularly important when administering the peptide orally where 

25 stomach acid can degrade the peptide. 

Another . embodiment of the present invention is to modify the 
amino acid sequence of a peptide to enhance its immunogenicity . 
This is done by modifying the natural peptide sequence to bind to 
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the Class II MHC receptor DR1 with superior binding affinity for a 
Class II MHC receptor DR1 than the natural peptide sequence. This 
modified peptide is considered a synthetic peptide. Alternatively, 
the sequence may be modified to have a greater inhibition of HA 
5 (3 06-318) binding to a Class II MHC receptor DR1. 

Many amino acid changes are acceptable in the formation of a 
synthetic peptide. The changes may be for similar types of amino 
acids such as leucine for isoleucine or they may be for diverse 
types such as tyrosine for lysine. 

10 

Materials and Methods: 

The structural homology model for the DR1 Class II MHC was 
constructed using the QUANTA molecular modeling package (vision 3.2, 
Molecular Simulations, Inc., Burlington, MA) with the CHARMM and 

15 Protein Design modules . After alignment of the sequences as 

described below, gaps and loops were energy minimized using 100 
steps of steepest descents minimization followed by 100 steps of 
adopted basis set Newton -Raphe son (ABNR) minimization. Large gaps 
were closed using a fragment database from a selected set of 

20 high- resolution crystal structures. The resulting structure was 

minimized in vacuo using 1000 steps of steepest descents followed by 
an additional 1000 steps of ABNR minimization. A distance related 
electrostatic function was used in all calculations with a 
dielectric constant of 1.0. Non-bound parameter lists were updated 

25 every 20 steps with a cutoff distance of 15. OA. Non-bonded 

calculations were performed using a shifted potential function 
between 11. OA and 14. OA. An extended atom set was used with only 
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polar hydrogen atoms specifically placed. There were no explicit 
hydrogen bond energy calculations performed. 

All peptides were initially modeled using QUANTA in an extended 
chain conformation and subjected to 500 steps of ABNR minimization. 
The resulting structures remained essentially in extended chain 
conformations. Individual peptides were manually docked in several 
different orientations into the binding cleft region of the 
minimized DR1 structure. The resulting bimolecular complex was 
subjected to 5000 steps of steepest descents minimization with 
non-bonded interactions updated every five steps. After 
minimization, bound peptides remained essentially in extended chain 
conformations. The lowest energy complexes for each peptide were 
selected for further analysis . 

The selected peptide and DR1 complexes and the minimized DR1 
model were subjected to the following molecular dynamics regimen: 
300 steps of heating to 300°K, 600 steps of equilibration at 300 °K, 
and 1100 steps of production dynamics. During this simulation, the 
DR1 Ca atoms were constrained in their starting positions. All 
non-bonded interaction parameters were as stated for the 
minimization procedure. The lowest energy structure during the 
course of the production dynamics was selected and subjected to the 
5000 step minimization procedure described previously with the Ca 
restraints removed. The resulting structures were used for the 
binding energy calculations and for hydrogen bonding analysis. 

Hydrogen bonds were determined using the QUANTA default 
parameters. Maximum allowed distances were 2.5A between a hydrogen 
and the acceptor atom and 3 . 3A between the donor and acceptor atoms . 
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The minimum angle allowed between any set of atoms forming a 
hydrogen bond was 90° . 

Competitive Inhibition Binding Assay: 
5 HA peptide (the influenza hemagglutinin 307-319 T-cell epitope) 

was labeled with 125 I . The labeled HA peptides were then allowed to 
interact with purified DR1 molecules during incubation to allow 
formation of peptide/DRl complexes. After incubation, the 
peptide/DRl composition was exposed to a native gel for 

10 chromatographic separation or passed through a spun column to 

separate labeled peptide/DRl complex and free labelled peptide. 
When unlabeled peptides were added before incubation of labeled HA 
peptides and DR1, and if the unlabelled peptides had capacity for 
binding to DR1 simultaneous with 125 I-HA, there was a resultant 

15 decrease in radioactive signal associated with the DR1 . The extent 
of this decrease directly related to the binding capacity of the 
unlabeled unknown peptide . 

Structural Homology Model for the DR1 Molecule: 
20 The structural homology model was created, the reference 

molecule being the crystal structure of HLA-aw68. The HLA-aw68 
coordinates and subsequent sequence were obtained from the entry 
2HLA in the Brookhaven Protein Data Bank released January 15, 1991, 
which is incorporated herein by reference. The sequence for the DR1 
25 molecule was for the a x domain was reported by Klein and for the (i x 
domain, the study reported by Todd et al . (Nature 329, 599 (1987)). 

The sequence alignment is based on Brown et al . (Nature 332, 
845 (1988)). The complete alignment and numbering scheme for both 

9 
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are seen in Figure 1. The Class II, /3j and Class I a 2 domains 
regions were conserved with some variations at the ends where the 
two MHC's have different loop regions. The fourth B-strand in the or, 
domain of HLA-aw68 (residues 30-38) is disrupted in the DR1 model. 
5 Only three residues are in a j3- sheet conformation, probably due to 
the inserted glycine at position 28 before the strand and the 
large deletion in the loop region immediately after the strand. The 
two alpha-helical regions are clearly maintained. Both helices have 
been -observed to be discontinuous in the Class I molecules and are 

10 similar in the DR1 model. The a x domain helix is long and curves 
from residues 49ar to 7€a without significant disruption. It is 
essentially a single continuous helix. However, the ct 2 helical 
region is broken into two separate helices as with the Class I 
molecules. A short helix (52-63) is separated from a longer helix 

15 (68-94) by a deformed region without secondary structure. This 

deformation is more pronounced in the DR1 model as opposed to the 
Class I molecules due to an insertion. 

Influenza Hemagglutinin Peptide with DR1 : 

20 The amino acid residues 307-319 of influenza hemagglutinin 

(Pro-Lys- Tyr-Val-Lys-Gln-Asn-Thr-Leu-Lys-Leu-Ala-Thr) make up a 
well -documented linear T-cell epitope which has been shown to be 
HLA-DR1 restricted. With the demonstration that the influenza 
hemagglutinin epitope (referred to as the HA peptide) binds DRl,it 

25 was chosen to be modeled into the binding cleft. 

The peptide was initially inserted into the cleft so that Leu 
11 HA was in the vicinity of the hydrophobic pocket. This allowed 
Asn 7 to be near the middle charged and polar groups of the cleft . 
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The remaining residue of the motif (Lys 2) was near the vicinity of 
the remaining charged and polar residues at the end of the cleft. 
The only adjustment to the starting conformation was a slight 
rearrangement of the terminal peptide proline and Tyr 3 to alleviate 
obvious bad contacts . 

After the energy minimization of the bimolecular complex, the 
total energy was reduced to 483 kcal/mol. This reduction in energy 
was accomplished by alleviation of several bad contacts and also be 
formation of several hydrogen bonds. The sticking feature of this 
mode is lack of hydrogen bonds in the carboxy terminal half of the 
peptide. Only one hydrogen bond is identified between the backbone 
carbonyl group of Leu 9 and the side chain of the j3, Asn 77. In 
contrast, the amino terminal half has eleven identified 
interactions. Four of these interaction involve the peptide 
backbone residues Tyr 3, Val 4, and Gin 6. The remainder involve 
the side chains of Lys 2, Tyr 3, Lys 5 and Gin 6. Interestingly, 
Lys 5 is involved in more interactions (three) than Lys 2 (only 2) . 
No interactions were observed as anticipated with Asn 7. Instead, 
it was the glutamine at position 6 donating a hydrogen bond to the a l 
Asn 62 . No interactions were observed for the amino and carboxy 
termini . 

HA-YK Peptide with DR1 : 

The binding of the HA-YK peptide (Ala-Ala-Tyr-Ala-Ala-Ala- 
Ala-Ala-Ala-Lys-Ala-Ala) to the DR1 model was tested . In aligning 
the peptide in the cleft, it was deemed logical to insert the 
tyrosine residue into the hydrophobic region of the binding cleft. 
The lysine would then be in position to interact with the 
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hydrophilic groups in the other half of the cleft. The resulting 
peptide orientation is the opposite that used for the HA and the CS3 
(defined below) peptides. With the peptide oriented as described, 
the final docking position for the peptide was unclear. The 
5 hydrophobic pocket is quite large, and, at least in this model, 

could accommodate the peptide tyrosine in a number of positions by 
sliding the peptide lengthwise through the cleft. However, 
repositioning the peptide also repositions the lysine. There were 
primarily two positions for the lysine: one with the lysine inside 
10 the cleft and the second with it outside. Of the two positions, the 
former was the lower in energy by 4 6 kcal/mol and had the greater 
number of interactions with the protein* (11 vs. 7). Thus, the 
preferred orientation of the peptide appears to be with the lysine 
inside the binding cleft region. 

15 

CS3 subunit Pilin Peptide with DR1 : 

The suspected T-cell epitope for CS3 pilus subunit 63-78 
(Ser-Lys-Asn-Gly-Thr-Val-Thr-Trp-Ala-His-Glu-Thr»Asn-Asn-Ser-Ala) 
was modeled with the DR1 molecule. The peptide was inserted with 

20 lysine inside the cleft in the hydrophilic region. This placed the 
Thr 5 in the center of the binding cleft and the tryptophane 
(residue 8) near the hydrophobic region. The resulting minimized 
model had ten interactions between the peptide and the protein, 
three interactions with the peptide backbone and five with the 

25 peptide side chains. The remaining two were with the amino terminal 
of the peptide. All of the interactions were in either the first 
three residues, His 10 or Glu 11 in the peptide. No interactions 



12 



WO 95/31997 



PCT7US94/05697 



were observed in the center of the cleft or residues four through 
nine . 

CFA/l with DR1: 

A peptide identified as CFA/l (colonization factor antigen) 
(Val-Gly-Lys-Asn- Ile-Thr-Val-Thr-Ala-Ser-Val-Asp-Pro) was prepared 
and an attempt was made to "fit" the molecule into the cleft of the 
DR1 . The lysine at position 3 prevented insertion of the peptide. 



10 



15 



Results : 

The peptides chosen to dock in the DR1 model are shown in Table 
1. The peptides were docked manually in several orientations into 
the DR1 model . The peptides were then tested in biological binding 
assays with the following results : 

Table I 



Peptide 


Molecular Model 
predicted binding 


Binding in the 
bioassay 


HA (influenza 
hemagglutinin) 


Yes 


Yes 


HA-YK (synthetic 
peptide) 


Yes 


Yes 


CS3 Pilin subunit 


Yes 


Yes 


CFA/l 


No 


No 



20 



25 



30 



Quantitative measurement of the inhibition of CS3 63-78 and HA 
306-318 as compared to controls is shown in Fig. 31. 

The binding energy was calculated as the difference between the 
final DR1 and peptide complex and the sum of the energies for the 
minimized DR and peptide models individually. The data is shown in 
Table II. 
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Table II. 



Peptide 


Protein 


Residues 


Sequence 


Binding 
Energy 
(kcal/mol ) 


HA 


Influenza 
hemagglutinin 


306-318 


PKYVKQNTLKLAT 


-283 


HA-YK 


synthetic 
peptide 




AAYAAAAAAKAA 


-216 


CS3 


CS3 pilin 
subunit 


63-78 


SKNGTVTWAHETNNSA 


-245 



CSSo; and CS6fe with DR1 

Colonization factor antigen IV (CFA/IV is an antigen on the 

10 surface of many enterotoxigenic E. coli one component of which is 
CS6 . CS6 has two major subunits and a number of minor subunits. 
Several peptides from CS6 have been sequenced and assayed for 
potential inhibition of radiolabeled HA (306-318 ) /DR1 complex as a 
measure of immunogenicity . The sequences of the subunits are shown 

15 in Table III. 



Table III, 



Peptide 


Amino Acid Residues 


Sequence 


CS6a6 


63-75 


DEYGLGRLVNTAD 


CS6o:7 


80-92 


I I YQ I VDE KGKKK 


CS6cy8 


111-123 


LNYTSGEKKISPG 


CS6S1 


3-15 


WQYKSLDVNVNIE 


CS6S2 


42-54 


QL YTVEMT I PAGV 


CS6£3 


112-124 


TSYTFSAIYTGGE 


CS6S4 


123-135 


GEYPNSGYSSGTY 


CS6E5 


133-145 


GTYAGHLTVS FYS 



These peptides were assayed for inhibition of radioactively 
labeled HA (306-318) /DR1 . The results are demonstrated in Fig. 32. 
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The foregoing description of the specific embodiments reveal 
the general nature of the invention so that others can, by applying 
current knowledge, readily modify and/or adapt for various 
applications such specific embodiments without departing from the 
5 generic concept, and, therefore, such adaptations and modifications 
should and are intended to be comprehended within the meaning and 
range of equivalents of the disclosed embodiments. It is to be 
understood that the phraseology or terminology employed herein is 
for the purpose of description and not of limitation. 
10 All references mentioned in this application are incorporated 

by reference . 
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We Claim: 

1. A method of preliminarily screening peptides for immunogenicity 
comprising the steps of: 

5 1) creating a molecular model of receptor DR1 Class II MHC and 

minimizing the model of the DR1 ; 

2) modeling a peptide to be tested and minimizing the model of 
the peptide; and 

3) testing fit of model obtained in step 2 into the model 

10 obtained in step 1 to produce a composite receptor/peptide model. 

2. A computerized model comprising a model of the DR1 molecule 
having fitted in a cleft therein a model of a peptide. 

15 3. A method of claim 1 wherein, additionally, the receptor/peptide 
model is subjected to computer-simulated heating. 

4. A method of claim 1 further comprising, assaying the peptide by 
competitive inhibition binding to a Class II MHC receptor DR1. 

20 

5. A minimized peptide capable of binding to a Class II MHC 
receptor DR1 and inhibiting the binding of HA (306-318). 

6. A synthetic peptide, wherein the amino acid sequence of the 
25 minimized peptide according to claim 5 has been modified to have a 

superior binding affinity for a Class II MHC receptor DR1 to form at 
least a portion of the synthetic peptide. 
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7. A synthetic peptide, wherein the amino acid sequence of the 
minimized peptide according to claim 5, has been modified to have 
greater inhibition of HA (306-318) binding to a Class II MHC 
receptor DR1 to form at least a portion of the synthetic peptide. 

5 

8. A synthetic peptide according to claim 6, wherein an amino acid 
has been modified from a charged amino acid to an uncharged amino 
acid. 

10 9. A synthetic peptide according to claim 7, wherein an amino acid 
has been modified from a charged amino acid to an uncharged amino 
acid. 

10. A synthetic peptide according to claim 8, wherein said 
15 uncharged amino acid is alanine. 

11. A synthetic peptide according to claim 9, wherein said 
uncharged amino acid is alanine. 

20 12 . A minimized peptide according to claim 5, wherein the sequence 
is selected from the group consisting of PKYVKQNTLKLAT , AAYAAAAAAKAA 
and SKNGTVTWAHETNNSA. 

13. A minimized peptide according to claim 5, wherein the sequence 
25 is contained in a CFA. 
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14. A minimized peptide according to claim 13, wherein the sequence 
is selected from the group consisting of DEYGLGRLVNTAD , 
IIYQIVDEKGKKK, LNYTSGEKKISPG , WQYKSLDVNVNIE , QLYTVEMT I P AGV , 
TSYTFSAI YTGGE , GEYPNSGYSSGTY and GTYAGHLTVSFYS . 



15. A vaccine comprising: 

a minimized peptide according to claim 5; and 
an immunologically acceptable carrier. 

16. A vaccine comprising: 

a synthetic peptide according to claim 6; and 
an immunologically acceptable carrier. 

17. A vaccine comprising: 

a synthetic peptide according to claim 7; and 
an immunologically acceptable carrier. 

18. A method of eliciting an immune response in an animal 
comprising administering said animal with the vaccine according to 
claim 15. 



19. A method of eliciting an immune response in an animal 
comprising administering said animal with the vaccine according to 
claim 16 . 

20. A method of eliciting an immune response in an animal 
comprising administering said animal with the vaccine according to 
claim 17. 
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